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describe the low-temperature supersolid "droplet crystal" phase, predicted for 
a broad class of soft-core interactions. We identify the conditions on the inter- 
particle interaction that render such intriguing phase possible, and outline propos- 
als for its observation. We argue this to be the prototypical supersolid, at least in 
the context of assemblies of ultracold atoms. 
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1 Introduction 

At the time of this writing, agreement has yet to be reached over whether the 
Non-Classical Rotational Inertia (NCRI) in solid 4 He at low temperature, first re- 
ported by Kim and Chan 1 , marks the first actual observation of the long sought 
supersolid phase of matter. The sheer weakness of the superfluid signal, the pre- 
dicted resilience of superfluidity (SF) in the metastable, overpressurized liquid^, 
the ostensibly important role played by extended defects, as well as some puzzling 
experimental observations (such as the concomitant occurrence of the NCRI with 
stiffening of the crystal™^, have so far precluded a straightforward interpretation 
of the phenomenology in terms of any conventional framework of supersolidity. 
First principle numerical simulations of crystals of 4 He, have yielded rather strong 
evidence that superfluidity does not take place in a perfect crystal™™", and cast 
serious doubts on microscopic scenarios based on Bose-Einstein Condensation 
(BEC) of point defects, such as vacancies or interstitials^l 
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While there exist dissenting views'^, it seems fair to state that consensus is 
growing that the ground state of solid 4 He is commensurate, i.e., free of point 
defects, and non-superfluid. Indeed, mounting experimental and theoretical evi- 
dence suggests that, if supersolidity occurs at all in helium, its character is likely 
to deviate significantly from the physical picture laid out in the seminal works of 
Andreev and Lifshits^, Chester^ and LegeetlP^. In particular, extended defects 
such as grain boundaries 12 or dislocations 13 might play a significant role in the 
phenomenology^ 4 -! 

In order to make progress toward a more complete understanding of the physics 
of solid helium, and more generally of the supersolid phase, a useful first step 
seems to be the unequivocal identification of a simple physical system, whose 
supersolid character can be established unambiguously, by means of controlled 
measurements as well as reliable, first principle calculations. The investigation of 
such a system may shed light on yet poorly understood aspects of the supersolid 
phase, which could in turn help resolve the controversy over solid helium, as well 
as facilitate the detection of supersolid behavior in other systems. The question is, 
of course, which other physical system, besides helium, can be opined to feature 
a supersolid phase, in a region of its phase diagram accessible by current technol- 
ogy- 

Among simple atomic or molecular condensed matter systems, helium has al- 
ways been regarded as the ideal candidate to display supersolid behavior. Due to 
the very low reactivity at ordinary conditions of temperature and pressure, helium 
atoms can be to an excellent approximation considered elementary particles; for 
the most abundant isotope ( 4 He), these particles have spin zero, and thus obey 
Bose statistics. On account of its low atomic mass, as well as of the weakness 
of the inter-atomic potential, liquid helium escapes crystallization at low tempera- 
ture, under the pressure of its own vapor. Solidification occurs if moderate pressure 
is applied, but atomic excursions away from lattice sites remain significant in the 
crystalline phase. 

The closest condensed matter system to helium, that may enjoy similar proper- 
ties, is molecular hydrogen (H2), also an assembly of Bose particles. Actually, the 
mass of a parahydrogen molecule is half of that of a helium atom, which should 
lead to even greater quantum effects. However, the attractive well of the interac- 
tion between two such molecules is about three times deeper than that between 
two helium atoms. As a result, liquid hydrogen crystallizes at a relatively high (14 
K) temperature, significantly above that at which BEC and SF are expected to oc- 
cur. Albeit effects due to Bose statistics are detectable in the momentum distribu- 
tion of liquid parahydrogen near freezing^, in general the behavior of condensed 
molecular hydrogen is much closer to that of a classical system than to helium. 
Indeed, while some theoreticaP^-^ and experimental 18 evidence of superfluidity 
(and even possible supersolidity^) in small clusters of parahydrogen (comprising 
around twenty molecules) has been reported, agreement is virtua lly un animous 
that bulk condensed parahydrogen displays no superfiuid behaviour 20 ' 2 H 

As it turns out, the physics predicted for solid 4 He (or, parahydrogen), as it 
emerges from first-principle quantum simulations, is essentially a direct conse- 
quence of the strong repulsive core of the interatomic potential at short distance 
(less than ~ 2 A). It is the repulsive core that determines most of the thermody- 
namic properties of helium and other quantum solids and liquids, witness the fact 
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that a very simple model of Bose hard spheres reproduces rather accurately the 
phase diagram of condensed helium^. 

Computer simulation studies of classical crystals, based on the Lennard-Jones 
potential (also featuring a hard core), suggest that, much like in a highly quan- 
tal solid like helium 7 , a uniform gas of point defects (vacancies or interstitials) 
is unstable against phase separation 23 . This suggests that the basic physics of the 
crystals is dominated by the strong interaction among particles, and that quantum 
effects cannot really alter it qualitatively (unlike in the liquid phase). One is there- 
fore naturally led to pose the following questions: can a supersolid phase of matter 
be underlain, in a many-body system, by a pair-wise interaction among elemen- 
tary constituents featuring a "soft" repulsive core at short distances ? And, if that 
is the case, do such physical systems exist in nature, or are they at least artificially 
realizable ? 



2 A different scenario: modified interactions and ultracold gases 

If one wishes to retai n the simple picture of "elementary" particles interacting via 
a static pair potential^ 24 ! 25 [ spatially confined ultracold gases provide an alterna- 
tive, in many respects more promising option to the controlled observation of the 
supersolid phase. It is now seventeen years since the first successful observat ion of 
BEC in Rubidium gas, cooled down to temperature in the nanoKelvin range 26 27 . 
Impetuous scientific and technological progress has rendered the field of ultra- 
cold atoms the playground of choice, to which one can turn to address outstanding 
questions in condensed matter and many-body physics. 

Dilute assemblies of cold atoms represent almost a textbook many-body sys- 
tem, a) virtually free from the imperfections and "background noise" that often 
mask subtle physical effects in a sample of condensed matter, and b) affording 
a degree of control that no naturally occurring solid or liquid can match. Last 
but not least, there exist a number of techniques whereby the interaction between 
atoms or molecules can be altered, making it virtually an adjustable parameter. 
The simplest example of such techniques, is the so-called Feshbach resonanc^\ 
allowing one to vary the strength of the (short-ranged) interaction between two 
atoms or molecules, and even reverse its sign (i.e., turned from repulsive to at- 
tractive, or vice versa). But then, if novel artificial inter-particle potentials can be 
created, the field is wide open for the search for a specific type of interactions, 
or broad class thereof, for which supersolid (or other yet unexplored) phases may 
occur. 



2.1 Dipolar systems 

As mentioned above, the repulsive core of virtually any known interaction among 
atoms or molecules, is what prevents microscopic scenarios of supersolidity based 
on point defects from occurring, as such defects have the tendency to cluster to- 
gether. On the other hand, there exist theoretical predictions of supersolid behavior 
in two-dimensional systems wherein pair-wise interaction feature a slow decay at 
large inter-particle separations, e.g., ~ 1 /r", with n < 3. 

The physical picture is radically different, in this case, from that based on BEC of 
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point defects, and it is based on the peculiar behavior of a purely repulsive dipolar 
system in the proximity of crystallization. In an ordinary first-order phase transi- 
tion, there exists a finite range of density (p,p + dp) wherein no homogeneous 
phase is fhermodynamically stable. Rather, two phases of different density (i.e., 
p and p + dp), coexist, separated by a macroscopic interface. It can be shown, 
however, that in the presence of long-ranged interactions (such as the dipolar) the 
energy associated to such an interface contains a negative term, diverging loga- 
rithmically in the thermodynamic limitPH As a result, ordinary coexistence is en- 
ergetically disfavoured. On approaching the transition from the low-density (e.g., 
the liquid) phase, the system may lower its free energy by embedding sufficiently 
large solid domains (i.e., macroscopic "clusters") inside the liquid. In the low tem- 
perature limit, two effects are predicted: the transition of the liquid to a superfiuid, 
and the crystallization of solid bubbles into a lattice superstructure, resulting in a 
global supersolid phase (in fact, a whole set of different such phases^). 

A well-defined way of testing the above prediction in ultracold gases, makes 
use of atoms or molecules possessing a finite electric dipole moment. These par- 
ticles can be confined to quasi two dimensions, by means of an external harmonic 
potential in the direction perpendicular to the motion (the so-called "pancake" ge- 
ometry). Upon aligning all dipoles in the direction perpendicular to the plane, by 
means of a strong external electric field, one can study a system of Bosons inter- 
acting via a purely repulsive potential^ of the form 1/r 3 . 

There is a subtle aspect in the derivation of the above predictions, however. For, 
the negative contribution to the surface tension is proportional to (5p) 2 ln(/?p 1 ' 2 ), 
where R is the droplet size. If the jump in density 8p at the phase transition is 
small (and quantum Monte Carlo simulations! 29 ! 30 ! 3 suggest that this may indeed 
be the case), the characteristic value of R required to observe the above scenario 
may be exceedingly large, for practical purposes outside the reach of realistic ex- 
perimental or even numerical setups. 

The remainder of this article focuses on a different class of interactions, fea- 
turing a soft-core potential at short inter-particle distance, one which does not 
grow arbitrarily, but plateaus to a finite value. Below, we review theoretical pre- 
dictions, based on first principle numerical simulations, of a system of such soft- 
core bosons, in two dimensions. The purpose is twofold: on the one hand, such 
a system can be shown to support a supersolid phase, whose physics turns out to 
be quite intuitive, arguably simpler than the scenarios that have been thus far con- 
sidered for solid helium. Secondly, because this type of interaction appears to be 
realizable in assemblies of cold atoms, by means of the so-called Rydberg block- 
ade, cold atoms may provide an entirely new, likely more direct pathway to the 
observation in the laboratory of a supersolid phase. 



2.2 Soft-core bosons 

In order to understand how a pair-wise interaction featuring a soft repulsive core 
at short inter-particle separation can lead to the appearance of a supersolid phase, 
it is useful to elucidate first the nature of the classical crystalline phase that occurs 
in a system characterized by this kind of interaction. To this aim, we consider the 
simplest soft-core pair potential, i.e., one that is equal to some energy V > if 
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Fig. 1 One-dimensional system of particles interacting via a soft-sphere potential, equal to V > 
if the distance between two particles is less than their effective diameter a, and zero if it is 
greater. As the distance between nearest neighboring particles approaches a, arrangement shown 
in (b) is favored over that shown in (a). 

particles are less than a distance a apart, zero otherwise. 

Imagine lining up particles interacting via this potential in one dimension, as 
shown in Fig. [T] As the density of the system approaches the point where each 
particle barely starts touching its two nearest neighbors, as shown in Fig.[TJa), the 
potential energy per particle discontinuously goes from zero to V, in correspon- 
dence to a classical crystal with lattice constant a. However, the system can lower 
its energy by one half, by allowing each other particle to overlap in part with one 
of its neighbor (on the left, for example), i.e., effectively doubling the lattice con- 
stant of the crystal, as each unit cell contains now two particles. On increasing the 
density, the number of particles per unit cell increases, while the lattice constant 
remains unchanged. Crucial to this effect, of course, is the finite energy "cost" 
associated to particle overlap. 

The subject of multiple occupation crystals (or "cluster crystals", as they are 
often referred to), has been extensively investigated in classical materials science, 
where they arise in the context of colloids and macromolecules^- 4 l 35 l 36 l 37 l In par- 
ticular, even though the qualitative argument offered above made use of the simple 
barrier potential, a quantitative criterion has been formulated, that allows one to 
predict the average number of particles K in the unit cell of the cluster crystal, 
based on general features of the soft-core interactional. 

How does quantum mechanics alter this simple picture ? In two fundamental 
ways: a) owing to quantum derealization, particles can "hop" to adjacent clusters, 
and b) an effective attraction among identical particles sets in in the case of Bose 
statistics, as exchanges of identical particles results in a lower kinetic energy (this 
is, of course, the same physics that gives rise to Bose-Einstein Condensation). 

Consider for definiteness a system of identical Bose particles of mass m, and 
let their interaction be described by a pair-wise potential, only depending on the 
relative distance of two particles, for which an effective range a can be defined. If 
all distances are expressed in units of a, and £ = h 2 /ma 2 is the energy unit, the 
(dimensionless) many-body Hamiltonian can be expressed as follows: 



where r (/ is the distance between ith and j'th particles. The only requirement on 
v(r), is that v(r — > 0) = v ~ l,As> where r s is the (dimensionless) mean inter- 
particle distance. In other words, the energy barrier associated to the repulsive 

core of the interaction, should not be much greater than the characteristic kinetic 

J32I 

energy 12 ^. 
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Fig. 2 Color online. Density of particles for a system of soft-core bosons in two dimensions, 
at low temperature. The nominal interparticle distance is r s = 0.5 in all four cases shown. The 
strength of the interparticle potential is D=\ (top left panel), D = 5 (top right), D=10 (bottom 
left) and Z)=100 (bottom right). The number of particles in each simulated system is ,/V=144. All 
lengths are expressed in units of the soft-core radius. 

In order to render the discussion more quantitative, we now focus on the phase 
diagram of a specific model of soft-core Bosons. We assume the simplest form of 
the pair potential, namely the "barrier" one discussed above 

V(r)=D@(r-l) (2) 

where is the Heaviside's function and the potential strength D is expressed in 
units of £ . As it tuns out, this simple model embodies all the physics of inter- 
est that we wish to present here. We postpone therefore the discussion of more 
more realistic interactions, until after establishing the main physical results for the 
barrier potential. Henceforth, we shall assume that the system is confined in two 
dimensions. 

The thermodynamic properties of a two-dimensional system of spin-zero bosons, 
interacting via a soft-disk pair potential, can be investigated by means of computer 
simulations, which have the advantage of being essentially exact, uncertainties be- 
ing reducible in principle to arbitrary degree, given enough compu ter time. In par- 
ticular, computer simulations based on the Worm Algorithm'2SI22l have emerged in 
recent years as a powerful methodology, affording both numerical accuracy as well 
as physical insight. It is particularly well suited to investigate Bose superfiuids, in 
that it affords the simultaneous, direct computation of the superfiuid fraction, us- 
ing the well-known winding number estimator of Pollock and CeperleyESl as well 
as of the one-particle density matrix. 

The phase diagram of a two-dimensional system of bosons interacting via 
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Fig. 3 Color online. One-body density matrix n(r) for a system of soft disks, at different values 
of the barrier strength D (corresponding to three of the four shown in Fig.[2). The temperature is 
taken sufficiently low that no significant change is observed on lowering it even further, within 
the statistical uncertainties of the calculation, i.e., these are essentailly ground state results. The 
nominal interparticle distance is r s = 0.5 in all cases shown. All lengths are expressed in units 
of the soft-core radius. 



the potential ^ has been recently investigated by Saccani et al^, by means of 
computer simulations making use of the continuous-space Worm Algorithm. We 
briefly review the basic physical results here. 

Fig.|2]shows the density of particles in the low temperature limit (typically, re- 
sults become temperature-independent for T < £ , i.e., they are representative of 
ground state physics), computed for a two-dimensional system of soft disks. The 
nominal inter-particle distance is r s = 0.5 is all cases shown; on the other hand, 
the strength D of the interaction Q increases hundredfold from top left (D = 1 ) 
to bottom right (D = 100). For weak interactions, the system is in a gas phase, 
the density being uniform (top left panel of Fig. ELj, but the formation of a solid 
crystal of clusters of particles, is clearly seen for D>2. The number of particles 
K per cluster remains remarkably constant (around 7 particles per clusters) as D 
is raised from a value of 5 (top right panel in Fig. [2) to 10 (bottom left), to 100 
(bottom right). 

The spatial extension of clusters decreases as D >> 1, i.e., they become more 
and more compressed. In that limit, particles in the same cluster tend to "pile up", 
in order to limit their overlap with other particles, and its associated energy cost, 
to those particles in the same cluster. Particle hopping to adjacent lattice sites, is a 
process that involves tunnelling under a potential barrier ~ KD, while the distance 
between nearest-neighboring clusters is ~ \/K. In a range of values of D (see be- 
low), this leads at low T to a superfluid transition, as particle hopping establishes 
phase coherence throughout the system. 

This is quantitatively shown in Fig. [3] which displays low-temperature results 
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Fig. 4 Color online. One-body density matrix n{r) for a system of soft disks, computed at 
different temperatures (in units of e ), for a value of the barrier strength D=5. The nominal 
interparticle distance is r s = 0.5 in all cases shown. All lengths are expressed in units of the 
soft-core radius. 

for the circularly and translationally averaged one-body density matrix, defined as 

n{r) = 4^Y I dQ S 11 ^ <^V + rW)> (3) 

where ~f is the volume of the simulation cell, / dfl represents angular integration 
and (...) stands for thermal average. Three different physical regimes can be iden- 
tified. At low values of D, the system is a uniform gas (D = 1), and consistently 
n(r) displays a slow power-law decay. At the opposite end, namely for strong 
inter-particle interactions, the system is an insulating droplet crystal; thus, n(r) 
decays rapidly, falling to essentially unmeasurable values for r > r s , i.e., outside 
the nominal inter-particle separation. This is because tunnelling between adjacent 
particles is exponentially suppressed by the height of the effective potential barrier 
seen by an individual particles traveling across droplets. 

For intermediate couplings, (i.e., D ~ 5), the one-body density matrix mir- 
rors the dual character of the system. On the one hand, n(r) decays slowly at 
long distances, in conformity with the presence of a superfluid regime described 
by Kosterlitz-Thouless theory 42 . Concurrently, it displays marked oscillations, re- 
flecting the underlying crystalline arrangement of the clusters. These oscillations 
arise microscopically from tunnelling of particles across adjacent clusters, and 
give rise to satellite peaks in the Fourier transform of n(r), namely the momentum 
distribution, which in cold atom assemblies lends itself to experimental imaging 
by means of time-of-flight measurements^]. Thus, the supersolid character of the 
system could be experimentally ascertained in a rather straightforward way. The 
emergence of off-diagonal long-range order, coupled with the solidlike oscilla- 
tions, is clearly shown by the temperature dependence of the n(r), displayed in 
Fig. [4] The superfluid character of the droplet crystal phase can be established by 
direct computation of the superfluid density 41 . A typical result is shown in Fig.pl 
for two simulated systems comprising different numbers of particles. As in any nu- 
merical study based on computations performed for systems of finite size, results 
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Fig. 5 Color online. Superfluid fraction p s , computed as a function of the temperature (ex- 
pressed in units of e Q ) for simulated systems comprising N = 144 (circles) and N=576 (dia- 
monds) particles. The nominal inter-particle distance is r s = 0.5, whereas the value of the cou- 
pling D=5. Statistical errors are smaller than the sizes of the symbols. 



must be obtained for several systems, comprising different numbers of particles, 
in order to carry out extrapolation of the physical estimates to the thermodynamic 
limit. The results shown in Fig. [5] are consistent with a superfluid transition in the 
Kosterlitz-Thouless universality class^l. It is noteworthy that the superfluid frac- 
tion does not saturate to a value of 100% in the T — > limit, consistently with the 
prediction of Leggett 1 1 for a superfluid system breaking translational invariance. 

Interestingly, it is possible to ascertain that, even in the insulating (e.g., non- 
superfluid) droplet crystal phase, individual droplets are superfluid, as shown by 
the statistics of permutation exchange cycles^. However, the transition of the 
system to a globally superfluid phase can only be established numerically in a 
well-defined range of values of D, as tunnelling is exponentially depressed with 
increased D. Because the computational methodology adopted here is a finite- 
temperature one, we cannot make a definitive prediction of a quantum phase tran- 
sition at 7=0, driven by the parameter D, between a supersolid and an insulating 
droplet crystal, based on the results of the simulation alone; however, our results 
are certainly consistent with that conclusion. 

The physical behavior of the supersolid droplet crystal phase is certainly remi- 
niscent of that of a Josephson junction array. Futhermore, for a system of identical 
bosons, one is naturally led to think of this many-body system in terns of a lattice 
Bose Hubbard model (BHM), which is known to display a superfluid ground state 
at low temperature^. It is important to note, however, that in this case the lattice 
is self-assembled, i.e., not imposed mathematically through space discretization, 
or experimentally (e.g., by means of optical lattices). It is this crucial difference, 
that ultimately allows one to speak meaningfully of a supersolid phase in this 
context. 
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2.3 Rydberg blockade 

The physics illustrated above would reain of merely "academic" interest, were it 
not because a well-defined procedure exists to realize a soft-sphere type interac- 
tion in an assembly of ultracold atoms, by means of a procedure known as Rydberg 
blockade (RB). Such a physical mechanism was initially introduced as a device to 
manipulate quantum information stored in collective states of mesoscopic ensem- 
bles' 4 -^, but it also has been recently proposed as a way to engineer a novel type 
of interaction potential between cold atoms. Specifically, the modified interaction 
"flattens off, and remains essentially constant below some characteristic "cut-off 
distance cffi. 

The soft sphere model should capture the main features of the interaction po- 
tential arising in the RB regime. However, it is necessary to assess the robustness 
of the main physical predictions illustrated above, by studying systems character- 
ized by more realistic model interactions, closer to those which could be realized 
in practice. In particular, one need assess any possible dependence on the physics 
of interest on the detailed features of the potential near the cut-off distance, as well 
as any influence of a long-range tail of the interaction itself. 

As it turns out, the main physical results, chiefly the presence of a supersolid 
phase, can be shown to be largely independent on details of the interaction, beyond 
the presence of a soft core. Thus, the droplet crystal phase ought be observable un- 
der relatively broad experimental conditions. For example, in the original work in 
which the droplet crystal phase was first predicted^, by means of computer sim- 
ulations like the ones described above, it was found that the same physics takes 
place in a two-dimensional system of bosons interacting through the following 
potential: 

V(r) « -j-J-g (4) 
r 5 +a J 

which also has a soft core at short distance, but a smooth behavior near the charac- 
teristic cut-off distance a, as well as a dipole-like long-range tail. The potential Q 
provides a realistic description of the specific type of interaction between pair of 
cold atoms that could be realized in the laboratory, via theRB 49 . All of this con- 
stitutes "circumstantial" evidence that the presence of a "plateau" at sufficiently 
short distances, namely a range within which the potential energy of interaction is 
either constant or varies very slowly as two particles are brought closer together, 
is the crucial feature that enables the formation of the droplet crystal phase, which 
turns supersolid at low temperature. Indeed, recent studies of two-dimensional 
Bose systems in which particles interact via a repulsive Yukawa potential show 
that no droplet crystal (nor any other supersolid phases) exist in that system^. The 
Yukawa potential has an infinitely strong repulsive core in the r — > limit, which, 
despite its very slow growth, suffices to prevent the occurrence of the physics dis- 
cussed here. Indeed, the basic physics of solids made of Yukawa bosons is the 
same as that of solid 4 He. It is also worth mentioning that no supersolid phase 
has also been thus far observed in numerical studies of purely repulsive dipolar 
systems in two dimensions. 

It has been recently conjectured^, that a necessary condition for the occur- 
rence of a supersolid phase (presumably of the "droplet crystal" kind), may be 
that the Fourier transform of the two-body potential go negative (i.e., that the po- 
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tential become attractive) in some range of k. Such a condition is not sufficient , as 
the purely repulsive soft-core interaction, whose Fourier transform is oscillating 
and indeed goes negative for some values of k, does not allow for a supersolid 
phase if the repulsive potential barrier is too high^M 

Now, obviously any realistic interaction displays a hard core at sufficiently 
short distances, as electronic clouds of different atoms are prevented by Pauli 
exclusion principle from spatially overlapping. Therefore, the relevance of any 
soft-core model ultimately hinges on a substantial difference in scales between 
the characteristic radius of the region in which the repulsive part of the potential 
softens, and the (much smaller) radius of the hard core. Because the softening of 
the potential at short distances is an artifical feature, arising from the RB, it seems 
conceivable that one may be able to tune the parameters in order to realize such a 
condition. 

On the other hand, the behavior of the potential at large inter-particle separa- 
tion is irrelevant to the existence of the supersolid phase described above, which 
has been predicted in theoretical s tudies in which different power law decays of 
the long-range tail were considered'^^MSI. 



3 Conclusions 

After decades of intense investigation, aimed at identifying a supersolid phase in 
solid helium, it now appears as if such a phase might be detected and explored 
more easily in a different context, namely in cold atom assemblies. Besides pro- 
viding the experimenter with a remarkably "clean" and controllable many-body 
system, the existence of techniques capable of fashioning the desired inter-particle 
interactions pave the way to the observation of novel phases of matter, difficult to 
achieve or non-existent in ordinary condensed matter physics. 

In particular, a physical mchanism known as the Rydberg blockade can give 
rise to pair-wise potentials between cold atoms, featuring a relatively "flat" re- 
gion at short inter-particle separations. As shown above, this spefici property can 
underlie a supersolid phase, consisting of a self-assembled crystal of superfluid 
droplets. Tunnelling of particles among adjacent droplets can establish phase co- 
eherence and give rise to a supersolid phase, whose most direct experimental sig- 
nature are secondary peaks in the momentum distibution. 

Computer simulations give evidence that such phase should be observable un- 
der relatively broad conditions, as it is insensitive to the detailed behavior of the 
inter-particle potential at long distances. The same physical behavior is also ob- 
served with different potentials featuring a flat region at short distance. Hard core 
potentials, on the other hand, do not lead to a supersolid phase, even those whose 
growth at short distance is slow, like the Yukawa potential. 

To our knowledge, no condensed matter system (either experimentally known 
or even toy model) has been predicted to display the supersolid phase illustrated 
here. The remarkable intuitiveness and simplicity of the supersolid droplet crystal 
phase makes it not only a promising candidate for a direct, unambiguous obser- 
vation of the supersolid phenomenon, but also a potential convenient prototype 
supersolid, allowing one to explore aspects of this intriguing phase difficult to ac- 
cess experimentally in other systems, such as helium. This suggests that a system 
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of soft-core bosons may warrant further investigations, aimed at exploring addi- 
tional aspects of the supersolid phase, which could render its identification in other 
condensed matter systems. For instance, the study of the excitation spectrum of the 
supersolid phase, with particular attention to the presence of two separate modes, 
corresponding to the simultaneous breaking of translational and gauge symme- 
tries, seems a worthwhile undertaking. 

Also of interest would be the study of supersolid mixtures, which may display 
intriguing de-mixing properties, as well as the motion of impurities through a su- 
persolid. 
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